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MODE COMPETITION IN THE QUASI-OPTICAL GYROTRON 
I. INTRODUCTION 

The quasi-optical gyrotron 1-10 (QOG) is a promising millimeter and 

q 

submillimeter radiation source that has already demonstrated power 
levels < 150 kV with pulse durations up to 13 ps at electronic 
efficiencies < 14%. Potential applications, ranging from heating of 
fusion plasmas to short wavelength radars, require stable CV operation at 
high power levels within a narrow frequency band. Because of the strong 
fields developed inside the cavity, the high density of longitudinal 
resonator modes, and the relatively long time of operation, nonlinear 
interaction among many modes will inevitably occur.In fact, it is 
the nonlinear coupling among cavity modes* 1 that will determine the 
existence of a final steady state, its accessibility from the initial 
small signal phase, and whether it is a single or a multimode state. 

In short, the evolution of the cavity fields can be described as 
follows. The injected electron beam initially excites all the cavity 
eigenmodes characterized by a start-up current below the beam current. 
At small amplitude each mode grows exponentially in time with its linear 
growth rate, unaffected by the presence of other modes. Soon one or more 
modes reach finite amplitude and start indicting through the induced 
modifications in the distribution of the Kctron beam. When this 
happens a single mode may nonlinearly suppress all other unstable modes 
and eventually dominate. Equally well, it may destabilize modes that are 
linearly stable. In some cases the fastest linearly growing mode may be 
overtaken by another, linearly slower mode, causing mode switching. The 
final steady state, if one exists, may involve more than one large 
amplitude modes. The point to be made here is that the existence of a 
final steady state, be it single or multimode, is determined nonlinearly 
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and cannot be predicted from the linear behavior. 

Finding the regimes in parameter space associated with the desired 
type of final equilibrium proves too costly for direct numerical 
simulation, given that the performance of the QOG depends on four 
independent parameters. It should be stressed that the final equilibrium 
cannot always be obtained by examining the linear stability of every 
possible single mode at saturation. 1 ^* Even if a saturated single mode is 
found stable to small perturbations by other modes, there is no guarantee 
that the system will eventually evolve to this state. In general, more 
than one stable equilibrium exists for a given set of modes, however, 
only one of them is accessible from the appropriate initial conditions. 
In the case that the final steady state is dominated by a single mode, it 
is not necessarily the mode with the largest linear growth rate. 
Consequently, all participating modes must be treated on equal basis 
until (and if) one dominates. Mode selection and accessibility are 
inherently nonlinear processes. 

In this work an analytic model of the multimode dynamics is 
developed by expanding the nonlinear current in powers of the radiation 
amplitude. A set of equations is then obtained for the slow time scale 
evolution of the wave amplitudes and phases. The strength of the 
nonlinear interactions enters through coupling coefficients of known 
analytic dependence on the gyrotron parameters. Numerical simulations of 
particle trajectories are then replaced by a set of first order ordinary 
differential equations, one for each participating mode. Not only is the 
computation time reduced by orders of magnitude but, equally important, 
the final equilibria can be predicted analytically from the coupling 
coefficients. In this paper we will mainly discuss cases with up to 
three interacting modes, coupled up to fifth order in the wave amplitude. 


2 







The coupled equations obtained here can be applied to the description of 
12 

phase-locked operation of the QOG through the injection of a small 

external signal. They can also describe phase locking through the use of 
13 

a prebunching resonator, via certain modifications resulting from the 
use of a prebunched injected electron distribution in place of a uniform 
one. 
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II. GENERAL FORMALISM 


The configuration for the QOG is shown in Fig. 1. The cavity 
fields are expressed through the vector potential A(r,t), 


where A is a superposition of eigenmodes of the Fabry-Perot resonator, 

1 - -i« t 

A(r,t) - - £ A m (t) U m (r) e m e y + cc , (2a) 

u " <r) ‘ YY r> c ° sk ™ x • <2b> 

In the orthogonal coordinates (x,y,z) with x along the resonator axis and 

z along the electron beam axis, the transverse profiles u (r) are given 

PQ 

by the Hermite-Gaussian functions 


W r) * Y <y) Y (z) ■ Y (z> * Y($>) cxp (-^) 1 


(3) 


polynomials. The wavenumber k m is equal to w „/c, where the cavity 
frequencies w take the values 

nc r 2( Pm + q m + 1} - 1 f m nc 

“m " “ [ m " 1 + --- tan W] = 7" (m + P m + «*.> ’ (4) 


where L is the cavity length and p m and q m are integers; the notation p m , 
q m is used for labeling frequencies, « m * w mpq’ Acc P r ding to the weak 
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coupling approximation, the slowly varying amplitudes A ffl (t) are complex 
quantities with magnitude A ffl (t) and phase t), 

i$_(t) 

= V*> e (5) 

The superposition (2), involving various frequencies « m , may also contain 
various modal structures p m q m for a given frequency (i.e. degeneracy). 


(a) Particle Equations 

Ve use the guiding center description for the transverse particle 
coordinates 

x * x„ + p cos 0 , y = y + P sin 0 , 

g g 

( 6 ) 

p * * p g* - p i si " 0 ’ p y ' p gy * P 1 cos 6 ’ 
where (x^.y^) and ^ p gx ,p gy^ denote the guiding center position and 
momentum, p^ is the magnitude of the transverse momentum, 0 is the 
gyroangle, p = Pj/8, is the Larmor radius, £? c - Q q /y is the relativistic 

cyclotron frequency, where Q q = |e|B o /mc, and the relativistic factor y = 
2 2 1/2 

11 + (p^/mc) + (p z /mc) J . By averaging the exact Lorentz force 
equations in the vector potential representation over the fast (cyclotron) 
time scale, the slow-time-scale nonlinear relativistic equations of motion 
are cast in the form 

du. 

— ■ E "n a n J i < 1 V> sin <*n 4 V cos Vg > (7a > 

at n 


d0 

dt 



r 


« J.(k p) 

Z n=I a "^V" C ° S 


(7b) 
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In Eqs. (7) time has been normalized to « o “* « (B q /y o ) *, length to 

k ■ c/» , a * |e|A /mc^ is the normalized radiation amplitude, u is the 
o o n 1 1 n 

normalized momentum u = p/mc * yv/c and the relativistic factor y = (1 + 
2 2 1/2 

u^ + u z ) . The angle is the relative phase between the m-th mode 

and the particle, = 0 - « m t, evolving in time as 


d* 


_ « Ji(kp) 

^ “i a " “v - cos 


( 8 ) 


where S« n = (« n - 2 Q /r o ) « Q _1 is the zero-order frequency detuning. The 
prime (') signifies the Bessel function derivative in respect to the 
argument. The evolution of y is found combining Eqs. (7a) to (7b), to 
obtain 


il a „ J i( k n ,) sin <*„ . V cos k n x g . (9) 

Two approximations are implicit in Eqs. (7)—(9). The guiding center drift 
has been ignored and the axial relativistic momentum u z is taken as 

constant. They are both justified for a n « 1, kV » 1, since dx /dt ~ 

2 

dy^/dt “ 0(a n ) and du z /dt " 0(a n /kW), which are much smaller than du/dt ~ 
d0/dt " 0(a n ). 

(b) Field Equations 

Substituting Eq. (2) into Maxwell's equations and ignoring second 
2 “ 2 2 ~ 

order derivatives 3 A n /3t « w n A n » tfie slow-time evolution of the vector 
potentials is given by 


6 






( 10 ) 


The second term in the right-hand side of Eq. (10) is zero by definition 
for the vacuum eigenmodes defined in Eq. (2b), thus, all the information 
for the slow evolution of the fields is carried in the nonlinear 
transverse current j^. Equation (10) can also be derived from the 
general expression governing coherent wave interaction, 1 * 



in the limit of standing waves inside the resonator (zero group velocity, 

2 2 

v g * 0) and near vacuum dielectric constant c (weak beam limit: /» o « 

1), (dc/d») n = 2» n . 

Multiplying both sides of Eq. (10) by exp(i&> n t) and taking the fast 
time average over t q ~ 2n/« o isolates the n th mode on the left-hand side. 
Introducing the normalized variables used in Eqs. (7), and letting a n * 
|e|A n /mc = a n exp (i$ n )» we have 


3a„ i 4n el 


U (r) 


3t 2« me 


V< r 5“„> 


273 J NL^ ; “n 


1 -o ito t 

m - J 0 dt W 1 ' 0 e 


(12a) 


(12b) 


Multiplying both sides by the mode profile U n and integrating over the 
cavity volume lead to 
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i 2n|e| 

2 3 

«/mc J V„ 
o n 


I i3r °„< r > %<"“„> • 


(13) 


3 2 

Above V R = J d r U n is the weighted cavity volume, and the 
phenomenological term « n /Q n was added to account for resonator losses due 
to boundary effects, where the cavity Q n is defined by 



(14) 


The current is expressed in terms of the nonlinear distribution 
function f^, 


j NL ( r .t) - -|e| n b (r) J^d9 ^ co S 9 f NL (r,u,t) . (15) 

Letting f NL (r,u,t) = n Q n^r^) f^u^u^zjt), substituting (15) in Eq. 
(13), and expressing r ± * (x,y) and u ± « (u x ,u y ) in the guiding center 
representation Eq. (6), the right-hand side of Eq. (15) becomes 


2n c |e| n 




‘ lu r L du l u l ~ cos0 J i <k n 0) £ NL (u i 


where again f NL (r»« n ) = ( 1/T 0 ) J 0 T dt f NL (r,t) ex P( iw n t >- Expression (16) 
is finally recast in terms of the injected beam current 
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h - I'I c( Vo J dx g J % "b ( VV • 


(17) 


(--^K--“<%<*>„>>• 


(18) 


The dimensionless, volume averaged current on the right-hand side of Eq. (18), 


y I 1 f T o itt t f 

- -7“ dte I dr U < r > W 1 ’ 0 

u z Vn T o J 0 J V 

- — — [ dr U ( 

*hV n Jv 


(19a) 


(r) W r; V 


is given, in terms of 


f„ L from Eq. (16), by 


<J NL<"„> > ■ 77 7 WV h p n (y g> “ s k n* 

bn z 


\j z h 1„ <2) J^ dS Io d “ 2 •fo d “ i “' L 7 “ se ‘m/W 9 ’*'V 


(19b) 


where I Q * lell^/mc^ is the dimensionless beam current (Budker parameter), 
and = /J dx^dy^ n^x^.y^) is the effective beam cross section. 
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III. NONLINEAR ELECTRON DISTRIBUTION 


The nonlinear electron distribution in phase space f NL (u 1 ,u z ,0,z,t) 
will be obtained in the form of an asymptotic expansion in powers of the 
small amplitudes a R « 1. Once the form of f^ in the multimode field (2) 
is known, the nonlinear transverse current j NL is obtained from Eq. (19). 
The evolution of f NL along the particle trajectories Eqs. (7) is given by 
the collisionless kinetic equation 


!fNL tu 

3f NL 

dO 

3f NL 

+ ^ 

3f NL 

du^ 

3t Z 

3z 

dt 

30 

dt 

3u z 

dt 


By formally expanding f^ and ordering in powers 

£ nl . £<°>, a £<».,J £<J> , 

one obtains a hierarchy of equations 


of a J , 


0 . ( 20 ) 

a ~ 0 (a n ) « 1 


L f ( j> = L. .Df^ 1 ) , 

O 1 

where 


L o ~ 3t + u z L + 2 c 30 1 


( 21 ) 


( 22 a) 


L 1 - £7 J l< k n-» 


e ‘*„ i_ 

3u. 


1 J,(kp) 


4 *n M 
aej 


k p 


( 22 b) 


Operator L q is the total time derivative L q » d/dt along the unperturbed 
orbits of electrons gyrating in the external magnetic field B q 
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e' - 0 - fi c T , u ± ' - u ± , z' = Z - —| * n ' * * n - S» n t ,(23) 

(the characteristic curves of L Q f = 0 in the theory of partial 
differential equations) where 0' ■ 0(t'), 0 ■ 0(t) and t = t - t'. 
Replacing the exact trajectories in the right-hand side of (21) by the 
unperturbed orbits (23), (integration along the characteristics), it 
follows that 

f^* = J^dT Df^ -1) [u'(t), 0'(t), z'(t) ] , (24) 

We exploit the fact that the interaction time of an electron inside the 
cavity, of the order of the transit time t^ = V/0 z , is much shorter than 
the characteristic growth time t c _1 = d(lna m )/dt for the fields. Keeping 
a n constant during the integration (24) yields 


f<j> 


(25) 



JjO^p) 1 3 

(k n p) 2u x 30 


f (j-D 


cc, 


where T j * ~ tjand we used 

V'j-P - - Vj-i ■ - Vj * *Vj • 

It is now convenient to change variables, introducing the set C Q , Cj, ••• 
, Cj, defined by 


z 



c r 






( 26 ) 
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where dCj - - (u 2 /yV) dx^ and ^j( T j c °) * ^j_i’ Successive iterations of 
Eq. (25) yield 


£<J). J^wv J o ... 

<- i)J 0 Jjwv jjwv JJ 



where f^ 


is the injected beam distribution, and is written as 



i(0 



iv C 


Jl(k n p) 3 
2 i 3u ± 


wi L) 

(k n p) 2u ± 30 J 


(28) 


cc . 


The quantity v r = (Vy/u z ) 6« n , the product of the frequency detuning times 
the transit time through the cavity is the most important parameter for 
the mode coupling coefficients. The integrations, Eqs. (27) and (28), are 
carried out in Appendix A. The expressions for f^\ j*l,2,3 are given by 
(A1)-(A3). For j > 4 the full expressions are getting rather complicated. 
However, one may exploit the existing ordering inside f^ in powers of 
the large parameter Y\ * (kW/y)(u^/u 2 ) » 1. It is then easier to find an 
approximation for f^ to any order in j, Eq. (A6), by keeping the 
dominant contribution. Additional simplifications occur during the 
gyroangle averaging in obtaining the nonlinear current, given by Eqs. (A7) 
and (A8). 
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IV. MODE COUPLING EQUATIONS 


Substituting expression (A7) for the nonlinear current inside the 
right hand side of Eq. (A8) and keeping terms up to fifth order in complex 
amplitudes a m = a m exp(i+ m ), the general mode coupling equations assume 
the form 


A 

E 

{i,j, 


zo b n 

( C m;n a m + C klm;n a k a l a m + C ijklm;n a i a j a k a l a m ) * 


(29) 


A mode of given frequency « n interacts through all possible frequency 
combinations satisfying the resonant condition 


^ij.-.mjn * <“i ± “j ± ••• ± w m > " “ n = 0 » (30) 

A 

denoted by the notation {i,j,...,m} implies summation over all 

permutations among i,j,...,m (except n) inside Eq. (29), for a given set 

of resonant frequencies (30). A negative index -m will imply -w m , -♦ m , 

-v and a * in place of «. * , \> and a„ respectively, 
m n mm m n 

The dominant contribution to the coupling coefficients is 

given by 


C 


( 1 ) 

m;n 


(?) “ 2 S «;n f”) ]> v »<*> J!«l V<1> 17- • 

zo 1 


(31a) 
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C kLU - (if 1,4 W “» [* ( Hl.in' * *kl«,„)] j’du,^ (-I -) 3 

® pZ 3 C 3 C> df 

- li* v z > f«i vv — C dc 2 w — C dc 3 \<v — ■ 

(31b) 


C ijklm;n * ft) ^ S ijklm;n exp [* ft^ijklmjn* + *ijklm;n)] J^du^— (jj—) 
r" f ® r^i ^ ® r^4 

J* v »<‘> f«i vv r - D 5 * W jr C-- iT W — ■ 


where 


v_(C) * J, (kp) h (C) e 


-iv C 


(31c) 

(32) 


*.. ■ (♦• ± ± ••• ± ♦ ) - ♦ is the combined slow phase, and the 

ij...m;n l j m n 

nonlinear filling factor is given by 


L 



2 1 


. 2 1 
£- 
2 


P- (> 



(33) 


The results appearing in Eqs. (29)-(33) apply for any spatial 
profile and velocity distribution of the electron beam, as well as any 
radiation modal structure in the open resonator. 
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V. THREE MODE INTERACTION 


Equations (29) are quite general, involving all possible resonant 
combinations of modes. In many cases it has been observed that only a 
small number of modes participates in the final equilibrium with finite 
amplitude. In this section the three mode synchronous interaction, 
satisfying the condition 

+ Vj * , or “k + “i a 2w m * (34) 

will be considered in more detail. This is a special case of the four 
mode interaction, with to m = u> n in Eq. (30). Hence the frequencies 
are symmetric around « m , \ = w m ~ dco, “i = “ m + dw » where dw is an y 
frequency interval allowed by the mode spacing in the cavity. The 
possible resonant frequency triads entering are 


W m ' “l - \ • 

- « m - «1 + «i , 

- “m - “k + “k ' 



(35) 


and all 3! permutations of each arrangement, while the possible 
arrangements for are given, setting i * k and j - 1, by 


- w m " “1 + “1 " \ * “k 

- w m - “k + “k " w m + w m 

- «_ - ~ 

m m m m m 

“m - “l - \ + “k " \ 

“m - “l - “k + W 1 - “l 


W m - “1 * “k + W m - W m 


(36) 
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with all the 5! permutations of each. 

There has been a tendency recently to employ the following 

renormalized parameters in gyrotron theory: amplitude F n " Wlo 3 " 
1 2 

E /B 3 , , detuning A = <2/3,^ )S« , interaction length u = 

n o oi n xo n 

k o W ( 3 lo 2/20 zo ) and current 1 - I 0 (8n/V o )(Q/Y o e io 3> where V o 85 (^2)k Q 3 LW 2 

and I Q = Ij ) (Ams) |e |/m e c 3 x 3 10^. Setting = « 2 » an( * 

separating real and imaginary parts, Eqs. (29) are recast in the form 


+ ZT F 1 - ReF l F 1 + F 1 F ® 2 ReG l<" + F 2 2 F 3 Re ( G 123 e " i# ) 
at m=l 

* l i h f " 2 Ret>n " * i F - 2 f 2 2 f 3 r '( d ” 123 •■**) ■ <37a> 


^ ^ F 2 - EeF 2 F 2 * 1 F 2 F m 2 EeG 2m * F 1 F 2 F 3 Ee ( G 231 e “) 

* l l ^ Fl2 ^ ReD21 " * ii F " 2 '» F 2 F 3 Re ( D »231 e ‘*) ' <37t> 


dF^ 3 a o ( 

— + —^ F^ « Rel^ F- + E ^ F ReG. ♦ F/ F, Re G. 
dt 2Q, 3 3 3 m=l 3 m lm 2 1 ' 1 


2 P „.r r e - 

312 e J 


£. 2, F 3 F l 2 F . 2 EeD 31. * F 2 2 F 1 Re ( D m312 •'**) - < 37c > 

1=1 m=l m=l 


where T is the complex linear growth rate. The evolution of the slow phase 
* = ^ + * 3 - 2^ is given by 
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a* , 

— - 52 + im r. + r, - 2 r 9 

dt k 1 i L ‘ 


- F l 2l 4ll* g 31- 2 G 2l) - F 3 2 i "( G 13* G 33 

F 1 F 3 

- 1, I, lm K. * D 31« - 2D 21m) F l 2 F » 

f .il F 3 F 2 2 f i*1 F 1 F 2 2 

' £1 " l Pj ^123* 


- 2 G 2 3 ) * F 2 2 l "( G 22 - 

• ”) * F 1 F 3 Ib ( G 231 
2 

ta (W **) - F 1 F 3 


° 12 " ° 32 ) 

e- “) 


«(W “) 


(37d) 


where SB * (tt^ + « 3 - 2«2)/« o . Note that only one combined slow phase 
appears for each triad. 

The dependence on the gyrotron parameters I, y and A is contained 
in the coupling coefficients. It will be shown now that when y » 
r o 3^ 0 /2, it can be eliminated as an independent parameter by proper 
scaling. The computation of I*, G, and D, carried out in Appendix B, 
yields 


r „ ■ Vl F <V • 

G nm - V3 < 2 «W • <38) 

B „„l ■ Vs * 4 °<V W • 

vhere t - <^> 2 <e i( , 4 / 2 S z0 2 >-y < Vo> and 2 s * 

I(0 io 2 /8Q)(k o W) 3 (0 lo /3 zo ) 3 [J 1 '(k o p o )] 2 . The nonlinear filling factors S y 
S^, and are given in (Bll). The quantities T, G, D .given for a cold 
beam by the integrals (B16)-(B18), depend on the desynchronism parameter 
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\> m alone, the.product of the interaction length u with the detuning 


putting the mode coupling equations in the final reduced form 


^ * F 1 - I { S l R ' r 1 F l * S 3 Sj F 1 F » 2 ReG 


lm + S 3 F 2^3 Re (®123 6 **) 


S 5 VlV R «°U. * s 5 Re (^ml23 ' < 41 “> 


I { S l Rer 2 F 2 ♦ s 3 i VS RcG 2n * S 3 F 1 F 2 F 3 Re ( G 123 * “) 


* S 5 f Z W T * Re ^21m * S 5 E/. 2f 1 F 2 F 3 Re K: 
1*1 m*l m«l 


, (41b) 


dF - 
— ♦ F, - I 

A .r J 


{ S l R ey 3 * S 3 R 3 F m 2 ReG 3« * S 3 F 2 2r 1 Re ( G 312 **) 


S 5 Z, ?, R 3 F 1 F m ReD 31m * S 5 E F . F 2 P 1 Re K,312 « J ' < 41 ‘> 


It follows, from (41), that the normalized start-up current is 1 = 
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3 ^ o the interaction length u, the 


The beam parameters I, 0 ZQ , 
cavity Q and the three frequencies combine into only 3 control parameters: 
I, Vj and y>2 (v^ is related to the other tvo frequencies through the 

resonant condition). The interaction length u does not appear explicity, 
but only implicitly through I and v. For simple density profiles, such as 
a pencil or a thin annular beam, the filling factors may also be expressed 
as powers of a single factor s, 

Sj * s j+1 , s * cos k Q x g h(y g ) , (42) 

with x and y the effective values of x and y . Thus, they could also 
g g g g 

be absorbed in the scaling factor, redefining 

V - s l , 

V = s 2 I . (43) 

Since I is given analytically in terms of the gyrotron parameters, it 
suffices to compute tables of T, G and D as functions of v m once, to 
completely know the coupling coefficients for any combination of the 
parameters I, y, P ZQ , 0^ and frequencies w m . Note that the beam spatial 
profile enters only in the nonlinear filling factor Sj; in case of 
complicated beam profiles Sj can be computed independently from Eq. (33). 
The discussion so far has been limited to a cold beam. In case of thermal 
effects the general expressions (A7) and (A8) are applicable and the 
control space will increase by two parameters, the pitch angle spread and 
the energy spread. 

Equations (37) and their reduced form, Eqs. (41), are the basic 
result in this paper, and describe three mode interaction up to fifth 
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order in magnitude in the QOG. The diagonal terms G nn describe the 

effects of nonlinear saturation (or self excitation) for each mode. The 

cross-coupling terms G nm describe a phase-independt. . nonlinear damping 

(positive or negative) among different modes. The interaction through the 

G , terms does involve the slow phase and describes mode locking effects, 
nmi 

The fifth order terms D nm ^ are necessary to stabilize the system in 
regimes vhere most of the third order interactions are positive and 
mutually destabilizing. They are also required to account for potential 
amplitude bistability, 1 ^ i.e., the existence of tvo equilibria of the same 
frequency but of different amplitude. 

Similar results, to third order in amplitude, have been obtained 
for the conventional gyrotron 11 , using a somewhat different approach. The 
present derivation of the coefficients is considerably simpler, since it 
requires a (j+l)-fold integration for the j-th order, compared to the 
(2j+l)-fold integration required in the single particle approach. 11 Mode 
coupling equations, identical in structure, first appeared in the 
treatment of LASER cavities 1 ^. 
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VI. RESULTS AND CONCLUSIONS 


The normalized linear growth rate ReT(v n ) is given in Fig. 2. 
Tables of the third order coupling coefficients G ( v n » v m ) f° r Gaussian 
profile modes and a cold pencil electron beam have been obtained by 
numerical evaluation of the expressions (B5)-(B7). The contour plots for 
the real and imaginary parts of G nm , and G 2 3 j within the regime 0 < 

v n , v m < 3 are shown in Figs. 3a t 3b, and 3c respectively (the 

coefficients ^^12 are given in terms of G i23» G 312^ V 1 ,V 2^ = G 123^ v 3’ v 2^ 

Dotted lines represent regions of negative values for G, corresponding to 

mutually stabilizing influence among interacting modes in that region. 

The mode coupling is destabilizing in the regions with solid lines, 

corresponding to positive G. Note the absence of symmetry in the coupling 

coefficients, where, in general, G t G . Contour plots for the fifth 
nm mn 

order coupling coefficients D ( v n » v m * v i) as functions of v m , for 
selected values of v n appear in Figs. 4(a)-4(d). 

According to the final number of participating equilibrium modes, 
Eqs. (41) demonstrate three types of equilibria: (a) single mode, where 
one mode dominates the other two (b) two mode equilibria, where one mode 
is of negligible amplitude and (c) equilibria among three modes of 
comparable amplitudes. The complete set of equilibria among a given set 
of frequencies is given by the zeros of the right-hand side of (42). 
Since F r are defined positive, only solutions with F^, F 2 » F^ > 0 are 

acceptable. It is well known that every stable equilibrium of a nonlinear 
dissipative system, such as the system represented by Eqs. (41), is 
associated with a region of initial conditions in amplitude space (basin 
of attraction) that eventually fall into this equilibrium. Which 
equilibrium the system will choose to settle in cannot be analytically 
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predicted, in the absence of invariants of the motion. This can only be 
done by integration of the set of ordinary differential equations (41). 

The simplest case, a single mode equilibrium for the n-th mode, is 
found from 

— = F n f ReT + ReG F n 2 + ReD nnn F 4 1 = 0 , (44a) 

n V. n nn n nnn n / 

— = ImT - F 2 f ImG + ImD_F 2 ) . (44b) 

^ t n n v nn nnn n ) 

Equations (44) are obtained by dropping the cross-coupling terms from 
(37); the combined phase 4 is irrelevant in that case. Acceptable, non¬ 
trivial solutions F n > 0 for the steady-state amplitude should always 

exist, since a single mode saturates through particle trapping. The fifth 

order saturation term D is therefore important on two counts. First, 
nnn 

it is necessary for saturation in regimes where the third order G Rn is 
positive and causes self-excitation instead of suppression. Second, it is 
required to account for amplitude bistability, when two acceptable 

equilibrium values F r j and F^ exist for a given mode, since the third 

1/2 

order nonzero solutions, F n = ± (IWG^) , can provide at most one 
acceptable solution F > 0. At steady-state the slow phase $ n varies at a 
constant rate, Eq. (44b); d+ R /dt expresses the nonlinear frequency shift 
from the linear frequency w n - 

Equations (37) are used to provide some examples of mode 

competition in a QOG driven by a cold pencil beam in a typical parameter 

regime, with beam current I. = 13 A, v = 1.146, 0, /0 = 1, V/X = 5, 

° b o lo zo 

cavity Q * 40,000 and cavity length L = 48cm, corresponding to a frequency 
separation d«/(i > 0 ■ 0.003. In Figs. 4(a)-4(c) we examine mode interaction 
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among the first three frequencies < « 2 < « 3 above w q , counting only the 
odd resonator modes o> n separated by du/o> o = 0.006. In Fig. 5(a) initially 
F 2o » F^ 0 » Fjq* The mode F 2 quickly reaches a quasi-steady state, which, 
however, is unstable. The mode F^, which initially grows much slowly than 
its uncoupled linear rate due to nonlinear suppression from F 2 , eventually 
overtakes and suppresses F 2 (mode switching in the cavity). In contrast, 
the steady state reached from F 1q » F 2q , F 3q , shown in Fig. 5(b), is 
stable since a large amplitude F^^ suppresses F 2 and F^ to negligible 
amplitude. However, the single mode equilibrium of F^ is not accessible 
from the usual start-up initial conditions, where all three modes have 
small comparable amplitudes F^ o = F 2q a « 1. Figure 5(c) shows that 
the final state in this case has F^ as the dominant mode. Note that F 2 
has been completely suppressed in Fig. 5(c), although its linear growth 
rate differs from F 2 by a few percent. Since F^ wins the competition with 
Fj and F 2 , we next examine the coupling of F^ with the next two 
frequencies F^ and F^. It is seen from Fig. 6 that the final steady state 
involves three modes with comparable amplitudes. Note that F 5 has been 
excited despite being linearly stable. Figure 7 shows the mode 
competition among the first three modes in a cavity of reduced length, L = 
24cm, and increased frequency separation dw/w Q = 0.012 among odd modes. 

In conclusion, the general analytic formalism for multimode 
interaction in the QOG has been developed in this paper. Applications of 
the coupling equations to interactions involving three main modes have 
demonstrated the ability to model nonlinear effects, such as nonlinear 
destabilization, suppression, and mode switching. A new scaling has been 
introduced, during the derivation of the coupling equations, that reduces 
the number of the control parameters. The ratio of operation to start-up 
current I, and the desynchronism parameters for the participating modes 
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become the only free parameters. The inclusion of coupling terms up to 
fifth order in amplitude was shown to stabilize the system in the hard 
excitation regime (where the third order coupling coefficients are 
positive.) The coupling coefficients up to fifth order have been 
tabulated for a cold pencil beam. These results will be generalized in a 
future paper to include thermal spreads in velocity, and arbitrary beam 
current profiles. The set of equations will also be expanded to include 
interactions among a large number of modes, by forming all the possible 
resonant frequency combinations. 
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Appendix A: Computation of the Nonlinear Current 


Substituting operator (28) inside (27) and performing the successive 
iterations yields 


kV . i(* + v O f C »f <0> 

Vfl.Ct) - - E a„ e j^X. n <W— 


(Al) 


v - v )c r c a r c i 

I - lX 


[**'■] 


V - V m>^ r c 1 ft 


) 3f<°> 

+ [ n » m ] j 3 Ujl 


(A2) 
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f (3) ( V e,t,t) . 0 3 E E I >„ a.* ®i*{ 

a r^i a r^2 


1( V V *i> i( V V V c 


+ £ -n, m, -1 j + £ -n, -m, 1 J + £ -n, -m, -1 ] + 

i<**i) i(v - v - v,K r c a r c i l 

[« 1 X. n -IjC 2 X. B -{jC 3 x 1 

+ [ -n, m, -1 j + [ -n, -m, 1 j + [ -n, -m, -1 ] + 


i( V V V i( V V V c 

e e 

- [ -n, -m, -1 ] 

i<V V V i<V V V< 

e e 

r i 1 8£<0> 


i r c i . a 

( c i r c i . l r^2 


(A3) 


where the explicit 0, t dependence is contained in %p n = 0 - « n t , 

y C,^) - •>; (^) «>* y g V« e 1V " C , xyc.u,) . y t.u/ 

zo 

X n <C,u l ) " X n (C,u l ) J 1 (~s7~) / J l(^T^) ("sr4 ’ X -n (C,u l ) * X n (C,u l ) 


0 o' 


In obtaining (A1)-(A3) we used the invariance: <J/ R ■*• v r C = const., along 
the unperturbed trajectories. The notation ( -n, m, -1 ] inside (A2)-(A3) 
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means an integral similar to the preceding fully written term except that 
a negative index -n in place of n implies -^ n , -v n> ^ a n * in place of »P n , 

v *-„• v etc - 

The dominant contribution in each integral comes from the 
derivative 3/3u^ acting on v m , 


3v « 

—- « — h(u,), vhere 

a.. 


k Vu. 2nW u, 

h(u.) = «-- » l . 

y u Xy u 


(A4) 


Therefore, the most important terms inside f^ are those with the 
(3/SUj^ derivative. Neglecting, for the same reason, terms like 3y/3u^ - 
3Uj^/3u^ - 1 « 3v m /3u^ one obtains from (A3) 


( Ul ,e,C,t) - £■) (^-) *! » n \ «!* 

zo 

*,> i < v n- v v* 3 r Ci 3 3f<0) 


+ £ -n, m, -1 J + [ -n, -m, 1 j + [ -n, -m, -1 j 

, X u n ivi t 

where v m <0 - Jj ^-J h m (C) e , 


(A5) 


One can show by induction, proceeding along the same lines, that 
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where { } implies all the permutations among the m^' s. 

When the j-th nonlinear piece of is substituted inside (15) 
and combined with cosG * (l/2)(e i0 + e _i0 ), the resulting expression 
contains even ©-harmonics, 1, e ±2i0 , ...» e ±U + l) i0 f f or j 0 dd f and odd 0- 
harmonics, e ±i0 , e ±3 * 0 , ...» e ±(j +1 ) i6 > f or j even . Consequently, coupling 
interactions involving even number of modes vanish completely during the 0- 
integration. The first nonlinear correction is of third order in amplitude, 
the second of fifth, and so on. It should be emphasized that this is a 
consequence of assuming a uniform in 0 initial distribution f^(u); 
however, for a prebunched distribution of the form f^ 0 \u,0) ■ f Q (u) + fj(u) 
e i0 + ... + cc, both even and odd terms survive and the nonlinear effects 
enter to second order in amplitude. Splitting the slowly varying ^ NL (« n ) 
component off the nonlinear distribution, by multiplying each f^ with 
exp(i« n t) and averaging over the fast time, substituting f NL (« n ) inside 
(19), and finally averaging over the gyroangle 0 and over the resonator 
volume, one obtains 


30 




<W“n» " - £ — “ 


ff b V n u zo t m i * m 2 * m 3 .V n} 


r<1) : . r(3) 


Vi ' a + C w/ a a a + ... + C VJ; a a . 

°* 1 5 n m j m l’ m 2’ m 3 ;n m l ra 2 m 3 m 1 ,m 2 ,m 3 ,...,m j ;n m 2 


C V^,n- -[‘K...,,.* * V-,- )] 

r" f a r c i 3 r^2 3 rS-1 df (0) 

b vz) b v -^ J.. dC2 v -^ > ^L- 


The quantities Aw m . v , # and S are 

m 1 m 2 ...m j ,n m. m 1 m 2 ...m j ,n nyy. .m jf n 

defined in a similar manner as in equations (30), (32) and (33) 

respectively. Of all the possible comoinations {m 1 ±m 2 ±m 3 ±...±m^] only 

those with (j-l)/2 positive and (j+l)/2 negative signs have survived the 

©-integration and contribute to for examnle. the term F-l.-m.-nl in 


©-integration and contribute to C VJ/ ; 
f^\ Eq. (A5), has vanished from c' 


for example, the term [-l,-m,-n] in 
(3) 

. . Also note that the complex 


conjugate terms in (A1)-(A6) drop out during the fast time averaging; they 
contribute to the j N i (-«_) component. 
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Appendix Bs Computation of the Three-Mode Coupling Coefficients 


The three-mode coupling coefficients in Eqs. (37) are computed from 
the general expressions in (31), keeping the frequency combinations (35) 
and (36), and yielding, 


r * K C ( *? n , G = K E C v m 
n -n;n nm -mm-n;n 

{-m,m,-n} 


r (3) 


K £ C —2—23;1 * G 312 = K E C -2-21;3’ G 231 = K E C 2-3-l; 


{- 2 ,- 2 , 3 } 


(- 2 ,- 2 , 1 } 


(2,-3,-1) 


w ml23 = K ^ ''-mm-2-23; 1 » 

{-m,m,-2,-2,3} 


r<5> 


= K E 

{-m,m,-1,1,-n) 


c (5) 

^-mm-ll-n;^ 


where K = 2ni 


(Bl) 


A cold pencil beam is employed for simplicity, where the velocity 
distribution f^ and the density profile n(Xg,y g ) are 8 ivGn by 


£ <0) - «<V « i0 ) ««, - u zo ) , n(x g ,y g ) - 5(x g -|) «y g ) . 
Znu lo 

The calculations are simplified by the following observation involving the 
derivatives of S(u^-u lo ). If A(u), B(u) are any functions of u satisfying 
d n A/du n « d n B/du n for every in, then 


J du A B ^ - J dug£ B 


i * dB I 

5 * - A sr l 


A(u 0 ) | , 
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and by induction 


r du A B £i * A(u ) fduB*. 
J rf„" 0 J Hu" 


jH / r \ i v ^ 

can see that -- I y, u. , J(k p) « -- e - h , where 

a „ n ' A n ' a .. n ° 


k V &ln 2 y 

>>1- 


(B2) 


Therefore, moving everything but the exponentials ex P( iv n O in front of the 
velocity integrals, one obtains 


r ‘ Sto li) U] .2 s. r, 
" 2 V„"zo \J [o1 1 " 


. V (J J ) S„ G 
o ) 'no mo' 3 nm 


no J lo J mo ) S 5 D nlm 
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(B5) 
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i,-l,l,-m,m) 

8 rS . 9 r c 4 . s . 

^L dt 4 V ”^L d ' 5V -"^ 8(ui ' uio) ’ (l 


iv t 

v n (C) - e n h (t) 


Ve limit the mode coupling among fundamental Gaussian modes where the 
transverse profiles Uqq are 

u 00 (y g’ z) = h 0 ( V h 0 (z) ’ h 0 (z) = exp (* ^2 ) * (B9) 

for all frequencies o» n . The normalized volume in that case is given by 


V n ' 2 / dx J dy I dz 


(‘ * h) 


exp(- = - k k 2 L W 2 , (BIO) 

' 9 n ° 


W(x) 


-c iv c 


For the pencil beam at x g = L/2, and according to Eq. (4), cos k n x g = 0 for 

n even and cos k x * 1 for n odd. In other words, the electron beam passes 
n g 

through a null of the electric field for modes with an even number of half 
wavelengths between the mirrors, or through a maximum, for an odd number of 
half wavelengths. The j-th order nonlinear filling factor Sj, Eq. (33), is 
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then reduced to 


r r 

S. « S „ m = 2 dx 
] V2 ...« r n J o gj__ 


% n( Xg ,y g ) (l * 


Jfi 


(j*Dy_ 


.(Bll) 


Since the frequency separation among the coupled frequencies » n is very 
small, it is important only when the detuning S« n , appears explicitly. 
Thus, one may set k R = k Q , letting ^'(^p) s * J 0 ’ V n = V = 


may notice that, according to (A4) and for » = « , 


d\> 


kW u. 


— = h(u.) = — — 


is independent of frequency. Since the S-function causes every derivative 
to be taken at u^ * u^ Q , one may employ v r , v^, as independent 
variables and use 


3u. ~ n o( 3v + 3v + 3v, ) * M 
1 n m 1 


to show that 
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where 
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It is clear from (B16)-(B18) that I\ G and D depend only on the initial 
values for the detuning parameters v nQ . Combining (B13)-(B15) with (B16)- 

(B18), converting the amplitudes a n into the normalized F n , a n = F n 

3 3 

(Y 0 @ i0 ), and absorbing the factors (r 0 B^ 0 ) in the coupling coefficients, 

finally yields 
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and, in general, for the j-th order coupling coefficient, 
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Defining the scaling factor 
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and the new normalized current 
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puts Eqs. (B19)-(B21) in the form given in Eq. (38). 
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Fig. 3 — Contour plots for the real parts of third order coupling coefficients (a) G.,, (b) G l23 and (c) G 2 
Dashed lines signify negative values. 
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(b) 

Fig. 4 (Continued) — Contour plots for the real parts of the fifth order coupling coefficients D„, m on the = 

const, plane for (a) v„ = 0.8 (b) = 1.9 (c) *>„ = 2.9. The contour plots in (d) are drawn for v m = . 


45 






Fig. 4 (Continued) — Contour plots for the real parts of the fifth order coupling coefficients D„ 
const, plane for (a) * 0.8 (b) *„ = 1.9 (c) •>„ = 2.9. The contour plots in (d) are drawn for y m 










(d) 

Fig. 4 (Continued) — Contour plots for the real parts of the fifth order coupling coefficients D.^ on the v n = 
const, plane for (a) v n = 0.8 (b) = 1.9 (c) = 2.9. The contour plots in (d) are drawn for = v,. 
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Fig. 5 (Continued) — Mode competition in the QOG for the first three modes above The time evolution of 
the amplitudes F in (a)-(c) corresponds to different initial conditions. 
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